Dynamical evolution of quintessence dark energy in collapsing dark matter halos 
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In this paper, we analyze the dynamical evolution of quintessence dark energy induced by the 
collapse of dark matter halos. Different from other previous studies, we develop a numerical strategy 
which allows us to calculate the dark energy evolution for the entire history of the spherical collapse 
of dark matter halos, without the need of separate treatments for linear, quasi-linear and nonlinear 
stages of the halo formation. It is found that the dark energy perturbations evolve with redshifts, and 
their specific behaviors depend on the quintessence potential as well as the collapsing process. The 
overall energy density perturbation is at the level of 10 -6 for cluster-sized halos. The perturbation 
amplitude decreases with the decrease of the halo mass. At a given redshift, the dark energy 
perturbation changes with the radius to the halo center, and can be either positive or negative 
depending on the contrast of dt<j>, dr<f> and <j> with respect to the background, where <j> is the 
quintessence field. For shells where the contrast of d r <j> is dominant, the dark energy perturbation 
is positive and can be as high as about 10 _J . 

PACS numbers: 98.80.-k, 95.36.-fx 
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I. INTRODUCTION 

The fast advance of different cosmological observations 
permits us to weigh the universe quantitatively [H, 0] • It 
is known that the universe is currently dominated by the 
dark energy component, which accounts for about 70% 
of the total energy density. Probing the nature of dark 
energy is thus one of the most challenging tasks in cos- 
mological studies. The simplest dark energy candidate 
is the cosmological constant with the equation of state 
w = — 1. While being consistent with observations, it 
runs into severe difficulties theoretically, namely the fine 
tuning problem and the coincidence problem. Various 
dynamical dark energy models have been proposed aim- 
ing to avoid such problems [1, 0, [f| . The allowed parame- 
ter space for dynamical dark energy models derived from 
current observational data is still relatively large @, 0] ■ 
Tighter constraints arc highly expected from future ob- 
servational data with much improved quality and quan- 
tity. To realize the goal, however, it is crucial to under- 
stand in detail the observational impacts from the dark 
energy component. Different dark energy models give rise 
to different expansion histories of the universe, which in 
turn influence differently on the formation and evolution 
of large-scale structures [1, H, EH ■ Besides, clustering 
can occur for dynamical dark energies themselves, and 
the clustering characteristics are highly model depen- 
dent. Therefore the observable imprints from dark en- 
ergy perturbations can potentially affect the constraints 
on the properties of dark energy. 

Dark energy perturbations associated with the linear 
evolution of the matter inhomogeneities have been stud- 
ied extensively [12|, |l3|. Considerable effects can occur. 
For instance, it has been found that with the proper in- 
clusion of dark energy perturbations, the constraints on 
the time-evolving equation of state of dark energy are 
weakened significantly due to their influence on the late 
ISW effects of CMB anisotropy Q1G3, OH. On the other 



hand, the dark energy clustering related to the nonlinear 
structure formation is less addressed. This is partially 
due to the numerical difficulties in solving the coupled 
perturbation equations for matter and dark energy in 
the nonlinear regime. However, many cosmological ob- 
servations, such as the abundance of clusters of galaxies 
and weak gravitational lensing effects on small scales, 
are directly related to the nonlinear structure formation 
[ItJ ■ It has been proposed that the recently reported cold 
spots seen on the CMB map may arise from the Rees- 
Sciama effect, the nonlinear version of the ISW effects 
[3 03 H3| • Therefore there is a need in understanding 
thoroughly the effects of dark energy on the formation of 
nonlinear structures. It is also of interest to investigate 
the corresponding clustering behavior of dark energy. 

There have been analyses on the formation of nonlin- 
ear structures, such as the virial radius and the average 
mass density of halos, taking into account the existence 
of the background dark energy [2l|, [H, [H, [24[ . Parame- 
terized dark energy perturbations and their correlations 
with the structure formation are investigated in some of 
the studies [22|, HE Hfl [27| • Recently, under the frame- 
work of spherical collapse of halos, the perturbations of 
quintessence fields induced by the structure formation 
are analyzed [28j, l29( . Different from the parameteriza- 
tion approach, in which the equation of state and the 
perturbation of dark energy are regarded as being in- 
dependent with each other, such analyses can present a 
consistent picture on the dark energy perturbation that 
often correlates with the evolution of the equation of state 
strongly. Dutta and Maor [28| mainly study the dark en- 
ergy perturbation in the linear regime of the halo forma- 
tion. Mota, Shaw & Silk [29[ consider the dark energy 
clustering in the process of nonlinear structure forma- 
tion employing a matched asymptotic expansion analysis 
method, in which, they adopt different parameter expan- 
sions for the linear, quasi -linear and nonlinear stages, re- 
spectively. In this paper, we present a numerical method 



2 



within the spherical collapse model that allows us to cal- 
culate the co-evolution of the halo formation and the dark 
energy clustering during the whole process, from linear 
to nonlinear epoches. In this approach, we do not need 
to treat the linear, quasi-linear and nonlinear evolutions 
separately. We analyze the quintessence dark energy per- 
turbation paying particular attention to its correlation 
with the background field evolution, i.e., the evolution of 
the equation of state of dark energy. 

The paper is organized as follows. In section 2, we 
present the model and the corresponding formulae. Sec- 
tion 3 explains in detail our numerical approach. Our 
results are shown in section 4. Section 5 contains discus- 
sions. 



II. THE MODEL AND FORMULAE 
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where the dot denotes the time derivative d/dt and the 
prime denotes the spatial derivative d I dr. 
From the Einstein field equations, we have 
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A. Two-components spherical collapse model 

In the framework of the spherical collapse of dark mat- 
ter halos, we consider the induced perturbation of dark 
energy. We start with the line element of a spherically 
symmetric system, which can be written as 



ds 2 
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The evolution of the quintessence scalar field, which is 
not coupled to the dark matter component except gravi- 
tationally, is governed by 
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--d a (V^9 a(3 d^) + ^- = 0, (6) 



where /i and R are functions of both time t and radius 
r. Note the corresponding R = a(t)r in the Fricdmann- 
Robertson- Walker (FRW) metric, where a(t) is the scale 
factor of the universe. The coordinate system is chosen 
such that it is comoving with the dark matter component. 
Thus the tt component of the energy-momentum tensor 
for dark matter is equal to the matter density, and the 
other components are zero. For the quintessence model, 
the energy-momentum tensor is 



t"P — 8 
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and 
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where V(</>) is the potential of the quintessence. The total 
energy-momentum T M „ can therefore be expressed as 



where -g = Det(g a p). 



B. Tolman approximation 

The perturbation of non-coupled quintessence dark en- 
ergy is expected to be small. Thus to a good approxi- 
mation, we can ignore the time-spatial term (/></)' in the 
right hand side of the equation (|5c[) . By time integrating 
its left hand side, we obtain a relation between /i and R, 
which follows 
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where / is a constant of integration and a function of r 
only, which reflects the local spatial curvature. Then we 
have 



Ttt = Pm + l<]> 2 + ^V 2 +V(<j>), (4a) 
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This is the well-known Tolman metric that is often uti- 
lized to compute the evolution of spherical objects. For 
a homogeneous universe, the Tolman metric reduces to 
the FRW metric and f (r) = —kr 2 . With the Tolman 
metric, the Einstein equations become 



R 1 / R 2 f 
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and 



D. The double exponential quintessence model 
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It is noted that keeping the cf>' 2 term formally does not 
affect our derivation of equations (|9aP and (|9b[) . On the 
other hand, to be consistent, we do not include it when 
numerically calculating the evolution of the metric R. 

The dynamical equation of quintessence ([5]) in the Tol- 
man metric is then written as 
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The equations (f9"aj) . (|9b|) and the equation (fit))) constitute 
the basic dynamical equations in our analyses. 



C. Spatial Curvature f(r) 

In the Tolman metric approximation, the term f(r) 
does not depend on time. Therefore we can calculate this 
quantity at any chosen time. In our consideration, we 
calculate f(r) at the initial redshift chosen to be z.- L = 30. 
Because the dark energy component is negligible at such 
a high redshift, the equations (|9a|) and (|9b|) can be further 
simplified. By combining the two equations, we obtain 



f(r) = R 2 \ Zi - 



2GM (r) 



R 



(11) 



where M (r) denotes the mass enclosed within the co- 
moving radius r, and 



dM (r) 
dr 



4nrp(zi,r) (R 2 R')\ 
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Therefore / (r) depends on the initial dark matter density 
profile. In this paper, we adopt the following density 
distribution 



p(ti,r) = p m (ti) 



2r^ 
3^2 



(13) 



where p m (U) is the initial background dark matter den- 
sity, S c is the initial central density contrast in the 
halo region, and a is related to the characteristic scale 
of the halo and will be discussed further later. For 
S c > 0/S c < 0, such a perturbation profile has an 
overdense/underdense zone in the central region and an 
underdense/overdense transitional middle region which 
compensates for the central region so that the overall 
perturbation is zero. 



In this paper, we consider the quintessence model with 
a commonly adopted double exponential potential 



V (0) = V (e Q ^ 
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where a and 3 are parameters representing the profiles 
of the two exponentials, respectively, and Vo reflects the 
height of the potential. In our analyses, we take a — 6 
and 3 = 0.1 as our fiducial model. We also vary a with 
a = 5, 6, 7 to see the effects of the potential on dark 
energy perturbations. The starting redshift is taken to 
be Zi = 30. The initial velocity of the <j> field is taken to 
be <j) = 0. The Vq parameter and the initial value of <f> 
are adjusted so that we have Qbe = 0.7 and u ~ —0.95 
at z — 0. Specifically, we first choose a pair of reasonable 
values of Vo and the initial </>, and evolve the background 
field to the present to calculate Qde and ui at z = 0. By 
comparing the obtained Qde and lu with CIde = 0.7 and 
lu ps —0.95, we adjust the values of Vq and the initial <j> to 
compute the evolution of the field again. Such a process is 
repeated until SIde — 0.7 and w w —0.95 are reached. It 
is known that the double exponential quintessence model 
has a late-time attractor solution, in which the final CIde 
and u) are mostly determined by the potential, i.e., the 
values of Vo, a and 8, but not very sensitive to the initial 
conditions of <j> [30]. Indeed, such a tracking behavior 
has been seen in our analyses. On the other hand, we 
start from a relatively low redshift (zt = 30), and choose 
4> = initially. Thus albeit not an extremely fine-tuning 
process, a certain level of adjustment on the initial value 
of (j) is needed in our analyses. 

Including the perturbations, we have 
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In Fig. [T] we show the evolution of the background 
equation of state w with V(f> = for a = 5, 6, and 7, 
respectively. The corresponding evolution of the time 
derivative of the background quintessence field, H<j), are 
shown in Fig. [2] It is seen that with the increase of a, 
the potential gets steeper and the H<fi term is larger. As 
a result, the equation of state u> increases. For a = 6 
and a — 7, uj > in the redshift ranges 5 < z < 14 and 
6 < z < 20, respectively. It will be seen later that the 
evolutionary behavior of the background field affects the 
dark energy perturbations significantly, thus in general, 
the two cannot be treated independently. 



4 




FIG. 1: The evolution of the background equation of state for 
the double exponential quintessence models. We take ft — 0.1, 
and the initial redshift zi = 30. The dotted, solid and dashed 
lines are for a — 5, 6 and 7, respectively. 




FIG. 2: The evolution of the background H(j>. The parameters 
are the same as in Fig. [1] 



III. NUMERICAL METHOD 



initial conditions for R are chosen to be R[zi) = a(zi)r 
and R — a(zi)r at Zi = 30, where a{zi) is the scale factor 
of the universe at Zj. 

With the obtained metric, we thus can analyze the spa- 
tial perturbation of dark energy. We rewrite the equa- 
tionfTU]) in the following form 
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R n 
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Equation (fTT)) can be regarded as a differential equa- 
tion of (j> with respect to r with the right hand side as 
the source of the spatial perturbations of cf>. We solve 
this equation iteratively at each time t. First, we insert 
the background 4>b g it) and 4>b g (t) into the right hand side 
of equation (fTT)) . With this source term, we solve 4>(r, t) 
at each r. Then for each obtained <f>{r, t), we find its 
corresponding <j)b g . This is done by looking through the 
evolution of the background field </>, and locating its 4>b g 
value at the place where <fi — (f>(r,t). Then at this time, 
different 4>bg value is obtained at different r. In the next 
iteration, this set of <fib g ( r ,t) are used in calculating the 
source term, which lead to a new set of (f>(r, t) . We repeat 
such an iteration process until the solutions on </>(r, t) 
converge. For each iteration, the shooting method from 
the center outward is applied to solve this spatial differ- 
ential equation (3l| . The boundary conditions are taken 
to be 4>{rb g ) = 4>b g and d(f>/dr\ r= o = 0. As shown in Eq. 
(|13p. we consider compensated matter density perturba- 
tions. Thus rb g is defined to be the shell inside which 
the average matter density perturbation is close to zero. 
The evolution of the shell at rb g follows that of the back- 
ground universe. 

In comparison with the method solving the partial dif- 
ferential equation directly, our numerical algorithm not 
only saves time but also avoids propagating numerical 
errors step by step in the time advancement. 



A. Iterative algorithm 



We first compute the local metric in the region of the 
spherical dark matter halo. We incorporate the local 
dark matter density and the dynamical evolution of the 
background quintessence field into our calculation. The 
dark energy perturbations are neglected in the metric 
calculations. Thus equation (f9"a[) can be written as follows 

5 + K")=- taG (^»- V *» ) ) (W) 

where bg stands for background. Equation (|9b[) has been 
used in deriving the expression of / given in part C of 
the previous section. 

Equation (|16p is an ordinary differential equation. We 
then solve it at different shells equally spaced in r. The 



B. Initial conditions 

In this study, we consider flat cosmologies with the pa- 
rameters being taken from WMAP5 results [32| . Specifi- 
cally, fl m = 0.249, h = 0.72 and er 8 = 0.787, where ft m is 
the dimensionless matter density of the universe, h is the 
Hubble constant in units of 100 km/s/Mpc, and as is the 
linearly extrapolated rms matter density fluctuation on 
top-hat scale of 8 Mpc/i -1 . The initial redshift is taken 
to be Zi =30. We consider overdense regions with typical 
scales of clusters and galaxies, respectively. For galaxies, 
we take M = 10 12 M©, and for clusters of galaxies, we 
take M — 10 14 Af Q . The corresponding linear comov- 
ing scales are calculated through M = (47r/3)il m p c rJ, 
where p c is the critical density of the universe at present. 
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FIG. 3: The local Hubble parameter R/R for different shells 
of dark matter. From upper to lower curves, r is decreasing. 



This gives rise to = 1.8 Mpc and r* = 8.4 Mpc, re- 
spectively. The initial perturbations over the considered 
scales are respectively set to be Si = 0.11 and Si — 0.044, 
estimated by employing the linear power spectrum of 
Eisenstein & Hu [33j ]. 

We now derive the relation between S c and a in equa- 
tion (|13|) and Si and r*. From the definition of Si, we 
have 



Si{r)AiTr 2 dr = 



where 



S,[r) = S,, - (1-^ 



(18) 



Then we obtain 



^ =Si 



In our analyses, we take S c = 5Si, therefore we have a = 
r*/1.27. The largest radius in our calculation is set to 
be r og = 6.5 Mpc for the galaxy scale and r^g = 30 Mpc 
for the cluster scale. Beyond the largest radius, all the 
relevant quantities merge into the background ones. 

For the curvature term f(r), as we discussed previ- 
ously, it can be calculated at Zi. Thus we have 



f = H 2 (z t )R 2 



2GM (r) 
Ri 



(19) 



We normalize R with the largest scales in our calcula- 
tions. 



C. Virialization 

In the general spherical collapse model, the virializa- 
tion process is not naturally included. It is known that 



TABLE I: Halo does not have a single collapse time in our 
model. This table shows A v i r of some shells at the time of 
their collapse. 



Zvir 


0.003 


1.07 


1.99 


3.57 


4.92 




148.96 


153.92 


162.33 


172.82 


179.21 



the violent realization and phase mixing play important 
roles in the final stage of the formation of dark matter 
halos. However, detailed understanding and modelling 
of the virialization process are still lacking. In [2^, [34[ , 
they fit the simulation results to describe the virialization 
stage. 

Here we adopt the following description for the virial- 
ization. The standard spherical model is used to calculate 
the metric evolution of R for each shell before its turn- 
around, denoted as R m ax- After the turn-around, the 
time derivative of R is modified such that it is approach- 
ing zero at the virialization radius R V i r . Specifically, we 
use 



R 1 / , R — Rmax 

~ = O I C0S '~E> 

RsCM " V Rmax 



Rv 



tt) + 1 , (20) 



where the SCM stands for the standard spherical col- 
lapse model. It is noted that in our modelling, the over- 
dense region has a r-dependent density profile. Conse- 
quently, different shells reach turn-around at different 
time. Thus Eq. (|20|) models the virialization for regions 
inside individual shells. In other words, the virialization 
starts from the inner region and gradually extends out- 
ward. To certain extents, we believe that our modelling 
grabs the basic picture of the virialization. 

Fig. [3] shows R/R for some shells, where we take 
Rvir(r) = Rmax(r)/l-8 instead of 2.0 from simula- 
tion results used in [2^, [35| . We also calculate the 
mean overdensity A„j r = p(r)/ P6 fl | z „ ir , where p(r) — 
M(r)/(4:TrR V ir(r) 3 /3) is the mean mass density within 
r and the mass M(r) contained inside r is computed by 
Eq. (I12p . and pb g is the background density at the viri- 
alization time. Here the virialization redshift z v ,- lr for a 
shell is defined to be the redshift when the shell collapses 
to the center in the case without considering the virializa- 
tion treatment. For different shells, the results are shown 
in Table |TJ It is seen that for regions where the virial- 
ization happens at high redshift, A V i r ~ 179, the well 
known value for fi m = 1. At later times, as the contri- 
bution from the dark energy component becomes signifi- 
cant, A v i r decreases. At z ~ 0, we have A« r ~ 150. Our 
results are consistent with those from the simple consid- 
eration of energy partitions [27j , indicating the meaning- 
fulness of our treatment for the virialization. 

In Fig. [U we show the time evolution of the relative 
density profile, where the horizontal axis is R/Rb g with 
Rb g being the R value at rj, ff and the verticle axis is the 
relative density p(R, z)j p og {z). As expected, the relative 
density gets higher and more centrally concentrated as 
time evolves. On the other hand, it is noted that our 



FIG. 4: The relative mass density profiles at different red- 
shifts. The solid, long dashed, dotted, dash-dotted and dot- 
dot-dashed lines are for the results at z = 20,5,2,1, and 0, 
respectively 

prescription of the virialization Eq. PU|) only accounts 
for the fact that R V ir / Rmax ~ 1-8 for a virialized shell. It 
does not include considerations about the detailed viri- 
alization processes, such as shell crossing, phase mixing 
and violent relaxation. Thus in our treatment, the evo- 
lution of R is basically an initial value problem, and the 
final relative density profile is closely related to its initial 
one. Since the NFW profile found from numerical simu- 
lations [13] must be at least partially related to the 
virialization processes, we do not expect that the NFW 
profile can be naturally obtained in our simple analyses 
presented here without adjusting the initial relative den- 
sity profile carefully. 



IV. RESULT 

In this section, we present the results on the calculated 
dark energy perturbations. 

Fig. [5] shows the time evolution of 5^/(1 + uj) for four 
spatial shells with r — 0, 0.25, 0.5 and 0.75, respectively, 
for the cluster-scale halo, where is the dark energy 
density perturbation, and r is in unit of the largest ra- 
dius (30.3 Mpc in this case) in our calculation. The cor- 
responding spatial variations of 8$/(l + uj) at different 
redshifts are shown in Fig. [5] For the purpose of easy 
comparison with the results of |29[, we follow them to 
show $0/ (1 + uj) instead of 5$. Under certain approxima- 
tions, it is shown in [29| that the quantity 6^/(1 +uj) can 
be directly related to the matter density perturbations. 
This motivates them to focus on 6^/(1 + uj). In our fol- 
lowing presentations, besides 5^/(1 + uj), we also show 
8$ at various places. We can see that the amplitude of 
<W(1 + UJ ) ls ~ 10 -5 — 10~ 6 except for the central region 
where 5^/(1 + uj) can reach as high as about 10 -4 . At 
1 + z > 20, the dark energy perturbations are positive 



FIG. 5: The logarithmic plot of the evolution of DE density 
contrast of four shells. The upper panel shows the positive 
part of 5^/(1 + u>) and the lower panel is for the negative 
part. The solid, dashed, dotted and dash-dotted lines are for 
r = 0, 0.25, 0.5 and 0.75 respectively. 




-2 1 1 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 6: The logarithmic plot of the spatial distribution of 
5^/(1 + uj) at of different redshifts. The upper panel is for 
the positive part of S^/(l + uj) and the lower panel is for 
the negative part. The solid, dashed, dotted and dash-dotted 
lines are for 1 + z — 22, 10, 1 and 0.12, respectively. 



everywhere in the calculated region, showing an initial 
concentration of dark energy in the region. With the 
evolution of the halo formation, the outer region appears 
as a dark energy void with negative 6^/(1 + uj), while 
the inner region (except for r = 0) has a positive dark 
energy perturbation. As time proceeds, the void region 
expands till very late time in the future when the dark en- 
ergy perturbation turns to be positive again everywhere 
in the region. 

In the following, we will analyze different contributions 
to the dark energy density perturbations in detail. Gen- 
erally, there are three terms to determine the total energy 
density of a scalar field, namely, the kinetic energy term 
~ 4> 2 , the potential energy term ~ V((p) and the term 
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related to the spatial variation of the field ~ 0/ 2 . For 
the very central shell with r = 0, we choose a boundary 
condition with cj>l = 0, and thus for this shell, only the 
kinetic and the potential terms exist. In Fig. the up- 
per panel shows the perturbations of the field <j> (dotted 
line) and its time derivative <p (dashed line) with respect 
to the background quantities. Because of its slower ex- 
pansion rate R than the background, the dragging term 
(R/R)4> is smaller in equation (JTUJ) . Therefore the evolu- 
tion of the field is faster, which leads to smaller value of 
<j) and the corresponding change of <f> comparing to those 
of the background quantities. The middle panel shows 
the corresponding energy perturbations. The lower panel 
presents 6$/(l+u), S%/(l+u) and (5$ + 8$)/(l+u) in 
logarithm scale, where the thick and thin lines correspond 
to positive and negative quantities, respectively. It is seen 
that the potential energy perturbation is always negative, 
while the kinetic energy perturbation can be positive or 
negative depending on the shape of the potential. At 
early times with z > 20, the kinetic energy perturbation 
is dominant, and thus we have 8$ > 0. At later times, 
the potential energy perturbation becomes the dominant 
term, which gives rise to the negative dark energy per- 
turbation, i.e., the dark energy void. It is noted that 
as the shell turns over, the term (R/R)<j> changes sign, 
and thus it is no longer the friction term, but a driving 
force for the local evolution of the field. It causes the de- 
viation of <f> from the background, and eventually leads 
to the domination of the kinetic energy perturbation and 
thus positive dark energy perturbations again at very late 
time with z < 0. It can be seen that the behavior of the 
dark energy perturbations is closely related to the form 
of the potential of the field, which in turn determines the 
evolution of equation of state of the background dark en- 
ergy. Fig. [8] shows the time evolution of the dark energy 
perturbations of the central shell for different potentials 
with a = 5 (dotted line), 6 (solid line) and 7 (dashed 
line), respectively. The dependence of 8$ on a is clearly 
seen. Therefore the dark energy perturbations and its 
background E.o.S are in general highly correlated, and 
cannot be treated independently. 

For shells between the center and the outer boundary 
the spatial variation of the field <f>l contributes signifi- 
cantly to the dark energy perturbations. Because of the 
homogeneity of the background, the spatial perturbation 
term always contributes positively to dark energy fluctu- 
ations. In Fig. [9j we show different parts of dark energy 
perturbations for the shell with r = 0.25. The upper 
panel shows the perturbations of <f> (dotted line) and <p 
(dashed line), respectively. In comparison with those of 
r = (the upper panel of Fig. [7]) , these perturbations are 
smaller. The middle panel shows the potential (dotted 
line), kinetic (dashed line), spatial (dash-dotted) terms, 
and the total (solid line) of dark energy perturbations. 
For this shell, the spatial term is the major perturba- 
tion term over a large redshift range with < z < 20. 
The lower panel shows the corresponding quantities over 
(1 + oj) in log-scale. It is seen that for < z < 20, the 
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FIG. 7: The evolution of the central density contrast. Upper 
panel: The evolution of A<f> (dotted) and A(j> (dashed); Middle 
panel: The evolution of 5f x 1CT 6 (dashed) and &\ x 10~ 6 
(dotted) and (Sf + Sjf) x 10~ 6 (solid); Lower panel: The 
evolution of logio\df /(l + a>)| (dashed) and logio\5^ / (1 + u)\ 
(dotted) and log l0 \(S$ + 5%)/{l+u)\ (solid). 



spatial term is much larger than the other two terms, and 
can be as large as about 10~ 5 . It is also noted that <j)l 
and the perturbation of <fi are correlated with each other. 

From Fig. [6l we see that the peak perturbation occurs 
at r ~ 0.1 at high redshifts (solid and dashed lines). It 
gradually moves outward, and eventually reaches r ~ 0.2. 
We find that the peak position is in good accordance with 
the shell within which the virialization happens. Such a 
shell corresponds to the position of maximum </>/. 

Fig. 1101 shows the average dark energy density pertur- 
bations within r — 0.25 (solid line), 0.5 (dashed line), 
0.75 (dotted line) and 1 (dash-dotted line), respectively. 
Firstly, we see the notable high peak at z ~ 10 for the 
solid curve, which is contributed by the spatial deriva- 
tive term </>/. As we average over larger regions outward, 
the dark energy void balances its concentration in the in- 
ner region and results smaller amplitude of dark energy 
perturbations. The inserted region of the plot shows the 
zoom-in pattern of the corresponding perturbations in 
very late times. The spatial variation of the dark energy 
perturbations can be relevant in considering the ISW ef- 
fect in a collapsing overdense region, for which we will 
investigate further in our future studies. 

To see the dependence of our results on the chosen 



8 




0.1 1 10 

1+z 



FIG. 8: The evolution of 84, x 10" 6 for r = 0. The dotted, 
solid and dashed lines are for a = 5, 6 and 7, respectively. 

initial density profile, we also perform the analyses with 
a compensated top-hat initial density distribution. The 
mass contained in the overdense region and the central 
density perturbation are set to be the same as those of 
the compensated Gaussian profile of Eq. (fT5)) . Thus in 
the top-hat profile, the overdense region is narrower than 
that of the Gaussian one. Fig. QT] shows the spatial dis- 
tribution of log 10 | <5<jf, / ( 1 +(jS)\ at z — with the solid and 
the dashed lines for the top-hat and the Gaussian pro- 
files, respectively. It is seen that the overall behaviors of 
the two are similar, with positive dark energy perturba- 
tions in the inner regions and negative ones in the outer 
regions. The narrower range of the positive perturba- 
tions in the top-hat case is closely related to its narrower 
matter overdense region than that of the Gaussian case. 
Since the profile of the initial matter density perturba- 
tions affects the detailed collapsing process, the specifics 
of dark energy perturbations are expected to depend on 
the initial matter density profile. However, the overall 
level of the dark energy perturbations is anticipated to 
be most sensitive to the total mass contained in the over- 
dense region. This can be seen from the following discus- 
sions. On the other hand, it is worthwhile to investigate 
in detail the dependence of the dark energy perturbations 
on the initial mass density profile, another direction that 
we can pursue in the future. 

We now compare dark energy perturbations induced by 
the formation of different sized halos. Fig. [T2l shows the 
dark energy perturbations at r — of the galaxy-sized 
halo with M = 10 12 M Q (solid line) and the cluster-sized 
halo with M = 10 14 M (dashed line). It is seen that 
the first transition from positive to negative dark energy 
perturbations occur at z ~ 20 for both cases. As we dis- 
cussed previously, this transition reflects the transition of 
the field from the fast evolutionary stage with a steep po- 
tential to the slow evolutionary stage with a relatively flat 
potential. Thus it is mainly determined by the behavior 
of the field potential itself, which explains the similarity 
seen in the two cases. On the other hand, the second 
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FIG. 9: Similar plot as in Fig. but for r = 0.25. Upper 
panel: The evolution of A<f> (dotted) and A(j> (dashed); Middle 
panel: The evolution of Sf x 1CT 6 (dashed) and 6^ x 10~ 6 
(dotted) and deformation energy density 5l x 10~ 6 (dash- 
dotted) and (5* + 5$ + 8%) x 1CT 6 (solid); Lower panel: The 
evolution of logio\Sf /(l + a>)| (dashed) and logio\d$/(l + u)\ 
(dotted) and logio\5$ / (I + lu)\ (dash-dotted) and log iq\(S$ + 
«£ + tf)/(l + w)| (solid). 
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FIG. 10: The average 84 x 10" 5 . The solid, dashed, dotted 
and dash-dotted lines are for the results averaging within r = 
0.25, 0.5, 0.75 and r = 1, respectively. 
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FIG. 11: The spatial distribution of logio\5,j, / (1 +w)| at z = 
0. The solid and dashed lines are for the compensated top- 
hat and the compensated Gaussian initial density profiles, 
respectively. 
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FIG. 12: The evolution of Zogi |^/(1 + w)| for r = 0. The 
solid line and dashed line are for the galactic halo and cluster 
halo, respectively. 



transition from negative to positive perturbations in late 
times is closely related to the formation of the halos. It 
happens earlier for galaxies than that of clusters. This 
is because of the higher density perturbations initially 
set for the galaxy-sized halo, which leads to its earlier 
formation than that of the cluster. For the overall ampli- 
tude of dark energy perturbations, it is about an order 
of magnitude larger for the cluster case than that of the 
galaxy case. Thus it is mainly the mass of the system 
rather than the initial amplitude of the dark matter den- 
sity perturbation that determines the overall amplitude 
of the quintessence dark energy perturbation. This is 
consistent with the expectation that quintessence dark 
energy perturbations show their effects mainly on very 
large scales. 



V. DISCUSSION AND CONCLUSION 

In this paper, we analyze the quintessence dark en- 
ergy perturbations induced by the spherical collapse of 
dark matter halos. We develop an iterative algorithm 
which allows us to solve the field perturbations consis- 
tently from linear to quasi-linear and to nonlinear stages 
of the halo formation. We find that the dark energy per- 
turbations depend sensitively on the field potential, and 
thus on the background equation of state. In the lin- 
ear regime of the halo formation, the slower expansion 
near the overdense region leads to the local evolution of 
the field to be ahead of the background field. When the 
quintessence field is in its fast evolutionary stage, \cj>\ is 
accelerating quickly. The perturbation of the kinetic dark 
energy term is positive, and dominates over the negative 
perturbation of the potential energy term. Adding in the 
positive perturbation from the spatial term cx <f>' 2 results 
a dark energy enhancement everywhere around an over- 
dense region at the very early stage of the evolution. As 
the field enters the region where the potential gets flatter, 
the |A</>| term is smaller, and the potential dark energy 
perturbation takes over. This gives rise to a negative 
perturbation, namely, dark energy void, in central and 
outer regions where the spatial variation term 0/ is small 
in accordance with our chosen boundary conditions. On 
the other hand, for regions in between, the (/>/ term is 
relatively large and there are strong dark energy concen- 
trations. In the late stage when the considered region coF 
lapses and virializes, the friction term (R! /R' + 2R/R)cf> 
changes sign to become a driving source for the field evo- 
lution. Thus the time evolution of <j> is accelerated al- 
though the potential is nearly flat. This leads to the 
domination of the positive kinetic dark energy perturba- 
tion again in very late stages. The overall dark energy 
perturbation is at the level ~ 10~ 5 for cluster-sized halos 
and ~ 10~ 6 for galactic halos. 

The work by [28( is the earliest one to analyze the 
dynamical dark energy evolution in the framework of 
spherical collapse of dark matter halos. They employ 
a linearized numerical approach, and thus only consider 
the linear stage of the halo formation. They find a scalar 
field void associate with the overdense region, and fur- 
ther argue that it should be deepened in the nonlinear 
stage. While our results from the linear stage of the 
halo formation are in qualitative agreement with theirs, 
some differences are worth mentioning. First, in their 
calculation of the dark energy perturbations, the term 
associated with the spatial variation of the field <pf is not 
included. It is seen from our analyses that in some spatial 
regions, this term can contribute significantly to positive 
dark energy perturbations to form a dark energy concen- 
tration rather than a void. Secondly, our results (e.g., 
Fig. [7|) show that the field void is indeed deepened ini- 
tially as expected, but then gradually becomes shallower. 
This behavior can be understood by noting that in the 
regions where the potential of the field is relatively flat, 
the velocity of the field \<p\ decreases. As a result, the 



10 



background field gradually catches up the field around 
the overdense region, and |A<^>| decreases. Therefore the 
depth that the field void can go depends sensitively on the 
dark energy potential, and thus the background equation 
of state. 

Our study is mostly related to the study in 29]. 
In their analyses, they adopt the method of matched 
asymptotic expansions (MAEs), in which different ex- 
pansion parameters are used at different regions, mainly 
HR << 1 in the interior region and S m << 1 in the ex- 
terior region. The two are required to match with each 
other in the intermediate region. With such approxi- 
mations, they derive a relation between the dark energy 
perturbations and the peculiar_ velocity Sv and the av- 
erage matter density contrast 5 m . Different treatments 
for the linear, quasi-linear and nonlinear evolution of the 
halo formation are then applied to obtain the information 
of Sv and S m , which finally give rise to the dark energy 
perturbations. On the other hand, we employ a full nu- 
merical treatment. The metric evolution is calculated in 
the Tolman approximation with the effect of dark energy 
perturbations neglected. We introduce an analytical de- 
scription to mimic the virialization of dark matter halos, 
which allows to carry out the metric calculation all the 
way from linear to nonlinear stages of the halo forma- 
tion consistently. With the obtained metric, an iterative 
algorithm is used to compute the dark energy pertur- 
bations. Given different approaches, our results are in 
good accordance with those of Mota et. al. [2!|. Both 
have seen the transition for the dark energy perturbation 
from positive to negative and to positive again. However, 
some quantitative differences exist. Considering the cen- 
tral region, the second transition from negative to pos- 
itive dark energy perturbations happens later from our 
results. We notice that in Mota et. al. 29] the metric 
calculation in nonlinear stages is done by ignoring the 
dark energy component completely. In our analyses, the 
dark energy perturbations are neglected, but the back- 
ground dark energy component is kept in computing the 
metric evolution. 

For a test, we have done a calculation on metric by 
dropping the dark energy component. The result on dark 
energy perturbations is shown in Fig. [13] for the central 
region (dashed line) together with the result with back- 
ground dark energy component included in metric calcu- 
lations (solid line). It is seen that the second transition 
indeed occurs earlier for the dashed line. This basically 
reflects that the existence of background dark energy af- 
fects the epoch of halo formation, which in turn affects 



the behavior of dark energy perturbations. We also find 
that while both catch the contribution of the <fi/ term to 
dark energy perturbations, we obtain a larger amplitude 
of the term than that of Mota et. al. This is closely 
related to the different initial density profiles as well as 
the different descriptions of the virialization process in 
the two studies. This, to certain extents, indicates that 
dark energy perturbations are sensitive to the details of 
their perturbing matter distributions. 
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FIG. 13: The influence of background dark energy on dark 
energy perturbation. The solid line and the dashed line are 
for the results with and without taking into account the back- 
ground dark energy in metric calculation. 

Our algorithm can be readily applied to study dark 
energy perturbations of different models. Further inves- 
tigations on their cosmological implications and seeking 
for possible observational effects of dark energy pertur- 
bations, e.g., the ISW effects around superclusters and 
clusters of galaxies, are highly desired and will be car- 
ried out as our future tasks. 
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